ELSEVIER 


Available  online  at  www.sciencedirect.com 


Journal  of  Power  Sources  139  (2005)  21-29 

www.elsevier.com /locate /jpowsour 


A  random  resistor  network  analysis  on  anodic  performance  enhancement 
of  solid  oxide  fuel  cells  by  penetrating  electrolyte  structures 

Dong  Hyup  Jeon,  Jin  Hyun  Nam,  Charn-Jung  Kim* 

School  of  Mechanical  and  Aerospace  Engineering,  Seoul  National  University,  Seoul  151-744,  South  Korea 

Received  16  May  2004;  received  in  revised  form  19  June  2004;  accepted  13  July  2004 
Available  online  25  September  2004 


Abstract 

The  anodic  performance  enhancement  of  solid  oxide  fuel  cells  (SOFCs)  by  introducing  penetrating  electrolyte  structures  was  investigated 
using  a  random  resistor  network  model  considering  the  transport  of  electrons  and  ions,  and  the  electrochemical  reaction  in  composite  anodes. 
The  composite  anode  was  modeled  as  a  mixture  of  ionic  and  electronic  particles,  randomly  distributed  at  simple  cubic  lattice  points.  The 
dependence  of  the  anodic  polarization  resistances  on  the  volume  fraction  of  the  electronic  phase,  the  thickness  of  the  anode,  and  the  insertion  of 
various  penetrating  electrolyte  structures  were  explored  to  obtain  design  criteria  for  best  performing  composite  anodes.  The  network  simulation 
showed  that  the  penetrating  electrolyte  structures  are  advantageous  over  flat  electrolytes  by  enabling  more  efficient  use  of  electrochemical 
reaction  sites,  and  thereby  reducing  the  polarization  resistances. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

SOFC  is  a  high  temperature  fuel  cell  that  directly  converts 
chemical  energy  of  gaseous  fuel  and  oxidant  to  electrical  en¬ 
ergy  by  electrochemical  reaction  [1] .  A  typical  SOFC  is  com¬ 
posed  of  a  composite  anode  made  of  nickel  (Ni)  and  yttria- 
stabilized  zirconia  (YSZ),  and  a  composite  cathode  made  of 
lanthanum  strontium  manganite  (LSM)  and  YSZ,  attached 
to  either  side  of  a  thin  YSZ  electrolyte.  The  performance  of 
SOFC  is,  therefore,  sum  of  the  performance  of  each  com¬ 
ponent  comprising  SOFC.  The  ionic  conductivity  of  YSZ  is 
around  0.1  S  cm'1  at  1000  °C.  The  low  ionic  conductivity  of 
the  electrolyte  is  among  main  factors  that  limit  the  perfor¬ 
mance  of  SOFC.  With  advance  in  material  processing,  thin¬ 
ner  electrolytes  become  possible  and  the  problem  due  to  high 
ohmic  loss  has  been  largely  overcome.  Currently,  electrolytes 
made  of  Y SZ  are  used  for  intermediate  temperature  operation 
of  SOFC  at  about  800  °C. 
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While  thinner  electrolytes  generally  perform  better  aside 
from  mechanical  problems,  that  is  not  true  with  electrodes. 
The  performances  of  the  anode  and  the  cathode  of  SOFC 
are  not  only  governed  by  the  transport  of  oxygen  ions  and 
electrons,  but  they  are  also  governed  by  the  electrochemical 
reaction  (also  by  the  transport  of  the  gas  molecules  of  fuel  and 
oxidant).  The  electrochemical  reaction  is  generally  accepted 
to  occur  at  three  phase  boundaries  (TPBs)  where  all  the  com¬ 
ponents  required  for  the  electrochemical  reaction  (oxygen 
ions,  electrons  and  gas  molecules)  co-exist  by  the  contact  of 
ionic  phase,  electronic  phase,  and  gas-diffusing  pores.  While 
reducing  the  thickness  of  the  electrode  is  favorable  for  the 
performance  by  decreasing  the  ohmic  resistance  and  by  in¬ 
creasing  the  rate  of  gaseous  diffusion,  it  also  reduces  TPBs 
where  the  electrochemical  reaction  occurs,  thereby  increas¬ 
ing  the  polarization  resistance. 

Many  studies  have  been  conducted  to  correlate  micro- 
structural  parameters  and  performance  of  the  composite 
electrodes  in  order  to  provide  fundamental  information  for 
optimization.  Using  either  resistor  network  models  based 
on  regular  lattices  or  random  packings  [2-4]  or  continuum 
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Nomenclature 


4 

current  generation  (A) 

4 

current  flowing  from  particle  i  to  particle  j  (A) 

Apb 

current  generation  per  TPB  (A/#  contact) 

/ 

length  (p,m) 

h 

ion  path  length  or  relative  distance  from  the 
electrolyte  (#  ionic  particle) 

ntpb 

number  of  TPBs  (#  contacts) 

R 

resistance  (£2) 

RpA 

area  specific  polarization  resistance  (m  Q  cm2) 

RtA 

area  specific  total  resistance  (m  Q  cm2) 

Greek  letters 

K 

conductivity  (S  cm"1) 

YAVb 

total  TPB  length  (p,m) 

y>tPb 

total  number  of  TPBs  (#  contact) 

Ga 

conductance  between  particle  i  and  particle  j 
(S) 

0el 

volume  fraction  of  electronic  phase 

potential  (V) 

Subscripts 

io 

ionic  phase 

el 

electronic  phase 

P 

polarization 

t 

total 

tpb 

three  phase  boundary 

Q 

ohmic 

Fig.  1.  Three  penetrating  electrolyte  structures  proposed  in  this  study.  Elec¬ 
trolytes  with  penetrating  block  structures  (a),  rib  structures  (b),  and  cross- 
linked  rib  structures  (c)  are  shown. 

electronic  phase,  oxygen  ion  migration  in  the  ionic  phase, 
and  the  reaction  at  TPBs  formed  by  the  contact  of  the  two 
phases  has  been  developed  to  quantitatively  assess  the  an¬ 
odic  performance  enhancement  by  the  penetrating  electrolyte 
structures. 


micro-models  based  on  the  statistics  of  random  packings 
[5-8],  criteria  for  the  volume  fraction  and  the  size  of  ionic 
and  electronic  particles,  and  thickness  of  the  composite  elec¬ 
trode  to  obtain  optimal  performance  have  been  suggested. 
Recently,  an  interesting  experiment  on  the  cathodic  perfor¬ 
mance  enhancement  by  depositing  relatively  large  YSZ  par¬ 
ticles  onto  the  surface  of  the  flat  electrolyte  has  been  reported 
by  Herbistritt  et  al.  [9].  After  a  numerical  analysis,  they  pro¬ 
posed  that  the  observed  better  performance  is  directly  related 
to  the  enlargement  of  the  surface  area  of  the  electrolyte  by 
the  deposited  YSZ  particles. 

By  extending  the  discussion  of  Herbistritt  et  al.  [9],  we 
proposed  penetrating  electrolyte  structures,  the  electrolyte 
structures  which  are  attached  to  the  flat  electrolyte  and  pen¬ 
etrate  into  the  electrode,  as  an  efficient  way  to  enhance  the 
anodic  performance.  Not  to  mention  the  increase  in  the  inter¬ 
facial  area  between  the  electrolyte  and  the  electrode  (primary 
reaction  sites)  by  the  penetrating  electrolyte  structures,  they 
are  believed  to  provide  less  resistive  paths  for  the  transport 
of  oxygen  ions  towards  the  interior  reaction  sites,  thereby  en¬ 
hancing  the  reaction  efficiency  of  TPBs  inside  the  electrode 
(secondary  reaction  sites).  In  this  study,  a  random  resistor 
network  model  considering  the  electron  conduction  in  the 


2.  Random  resistor  network  model 

Three  penetrating  electrolyte  structures,  proposed  in  this 
study,  are  shown  in  Fig.  l(a)-(c),  where  block  structures, 
rib  structures  and  cross-linked  rib  structures,  made  of  ionic 
phase,  are  attached  to  the  surface  of  the  flat  electrolytes.  In¬ 
troducing  the  bulk  electrolyte  structures  that  penetrate  the 
composite  anode  is  favorable  for  the  anodic  performance 
by  enhancing  the  oxygen  ion  transport,  thereby  facilitating 
the  efficient  utilization  of  TPBs  for  current  generation.  The 
penetrating  electrolyte  structures  are  also  favorable  for  the 
reduction  in  the  concentration  polarization  by  increasing  the 
limiting  current  density.  As  the  performance  enhancement  by 
the  proposed  penetrating  electrolyte  structures  is  purely  ge¬ 
ometrical,  such  geometrical  modification  can  be  superposed 
with  the  invention  of  new  materials  to  obtain  better  perfor¬ 
mance. 

From  the  mechanical  point  of  view,  the  penetrating  elec¬ 
trolyte  structures  are  believed  to  improve  the  reliability  of 
SOFC  significantly.  The  penetrating  electrolyte  structures 
reduce  the  volume  fraction  of  electronic  phase  in  the  com¬ 
posite  electrode,  thereby  relieving  the  problem  of  thermal  ex¬ 
pansion  mismatch  between  the  electrolyte  and  the  electrode. 
They  also  facilitate  interdigitated  close  contact  at  the  interface 
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Fig.  2.  A  random  resistor  network  model  for  a  composite  anode.  The  com¬ 
posite  anode  is  composed  of  equal-sized  ionic  and  electronic  particles,  placed 
on  a  flat  electrolyte  with  a  penetrating  electrolyte  structure. 


between  the  electrolyte  and  the  electrode,  which  significantly 
reduces  the  possibility  of  debonding  or  delamination.  The 
penetrating  cross-linked  rib  structures,  especially,  can  com¬ 
plement  small  buckling  resistances  of  the  electrolytes,  mak¬ 
ing  thinner  electrolytes  possible  [10]. 

The  composite  anode  of  SOFC  was  modeled  as  a  mixture 
of  the  ionic  particles  (white)  and  electronic  particles  (black) 
packed  between  the  electrolyte  and  the  current  collector,  as 
shown  in  Fig.  2.  Then,  the  resistor  network  was  constructed 
based  on  the  particle  network  of  equal- sized  spherical  parti¬ 
cles,  prescribed  randomly  as  either  the  ionic  particle  or  the 
electronic  particle  and  placed  at  the  lattice  points  of  a  simple 
cubic  lattice.  Though  the  resistor  networks  based  on  regular 
lattices  are  only  approximation  of  the  micro-structures  of  the 
real  composite  anodes,  such  an  approach  has  been  proven 
to  produce  reliable  prediction  when  appropriate  particle-to- 
particle  resistances  are  used  [2,4] . 

The  electrochemical  reaction  at  the  active  reaction  sites  or 
TPBs  requires  the  oxygen  ions  migrated  from  the  electrolyte 
towards  TPBs  through  the  network  of  the  ionic  particles,  and 
the  fuel  gas  molecules  diffused  from  the  channel  towards 
TPBs  through  pores.  Upon  completion  of  the  reaction,  the 
resultant  electrons  are  conducted  from  TPBs  towards  the  cur¬ 
rent  collector  through  the  network  of  the  electronic  parti¬ 
cles,  and  the  reactant  gas  molecules  are  transported  back  to 
the  channel  by  diffusion.  In  order  for  the  above-mentioned 
processes  to  happen,  the  ionic  particles  and  the  electronic 
particles  should  form  continuous  networks  (percolation)  as 
shown  in  Fig.  2,  assuming  the  continuous  pore  network  is  al¬ 
ways  available.  With  negligible  diffusion  overpotential  (con¬ 
centration  loss),  the  performance  of  the  composite  anode 
solely  depends  on  the  transport  of  ions  and  electrons,  and 
the  electrochemical  reaction.  In  this  case,  the  advantage  of 
the  penetrating  electrolyte  structures  in  shortening  the  path 
for  the  oxygen  ion  transport  can  be  easily  deduced  from 
Fig.  2. 

The  governing  equations  for  the  conservation  of  the  oxy¬ 
gen  ions  and  electrons  are  derived  from  the  Kirchhoff ’s  law 
of  current  [2] .  Current  conservations  for  a  particle  i  connected 


to  a  particle  j  are 

Ezio  \  A  Ao/,Ao  Jo\ 

Ii^j=  aiM  -<Pj) 

j=ionic  j=ionic 

=  ^  for  ionic  particle  i  (1) 

j=electronic 


£  '&/=  £ 

j=electronic  /=electronic 


=  —  ^  for  electronic  particle  i 


j= ionic 


(2) 


where  I^j  is  the  current  flowing  from  the  particle  i  to  the 
particle  j ,  Oij  is  the  bond  conductance,  (pi  and  cpj  are  the  po¬ 
tentials  of  the  particle  i  and  the  particle  j ,  respectively,  and 
is  the  charge  transfer  current  due  to  the  electrochemical 
reaction. 

The  charge  transfer  current  is  a  function  of  the  exchange 
current  density,  the  polarization  overpotential,  and  the  re¬ 
actant  and  the  product  concentrations  in  the  anode,  as  de¬ 
scribed  by  Butler- Volmer  equation.  At  low  current  density, 
the  Butler- Volmer  equation  can  be  approximated  by  a  linear 
charge  transfer  equation,  and  the  transfer  current  from  the 
electronic  particle  i  to  the  ionic  particle  j  is 


If 


lj-  =  aMx  -  vf), 

Act 


(3) 


where  Rct  is  the  charge  transfer  resistance,  ap  is  the  polariza¬ 
tion  conductance,  and  r)  is  the  polarization  overpotential,  r)  = 

<pel  ~(pio. 

Combining  Eqs.  (l)-(3)  results  in  a  single  form  of  the  cur¬ 
rent  conservation  equations,  written  for  a  particle  i  to  neigh¬ 
boring  particles  j  as 


V  =  V  -  <Pj)  =  °*  (4) 

j  j 


where  the  bond  conductance  crij  is  defined  as  <7i0_i0  for  the 
contact  of  two  ionic  particles,  aei-ei  for  the  contact  of  two 
electronic  particles,  and  ai0-ei  for  the  contact  of  an  ionic  par¬ 
ticle  and  an  electronic  particle. 

The  bond  conductance  is  dependent  on  the  geometry  of 
neck  formed  by  the  contact  of  two  spherical  particles.  Based 
on  Sunde  [2,3],  we  assumed  that  the  neck  perimeter  is  equal 
to  the  TPB  length  per  contact  of  dissimilar  particles.  The 
bond  conductances  between  the  same  particles  were  evalu¬ 
ated  using  the  neck  perimeter  /i0_i0  and  /ei— el ,  and  the  ionic 
conductivity  k\0  and  the  electronic  conductivity  /ce i,  as 

k'io^io— io  Kel/el— el  /c. 

^io— io  —  ~  5  ^el— el  —  “  •  w) 

The  bond  conductance  for  the  contact  of  dissimilar  particles 
is  determined  by  the  sum  of  ionic,  electronic,  and  polarization 
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resistances,  written  as 

/  i  i  i  v1 

CTio-el  —  (  - - 1  I-  - -  )  .  (6) 

\  2crei_ el  J 

where  the  polarization  conductance  ap  is  defined  by  the  po¬ 
larization  conductivity  /cp  times  TPB  length  /tpb,  as 

=  ^P^tpb-  (7) 


ence  of  unity  between  the  top  and  the  bottom  boundaries  (cpt op 
—  (pbot  =1),  the  distribution  of  the  potential  in  the  composite 
anode  was  determined.  For  the  solution  of  the  system  of  gov¬ 
erning  equations,  the  method  of  successive  substitution  was 
used  with  Gauss-Seidel  iterative  solver.  The  convergence  of 
the  potential  was  judged  by  the  divergence  of  the  current  for 
all  percolated  particles  i  to  be  smaller  than  a  predetermined 
value  as 


For  the  consideration  of  the  penetrating  electrolyte  structures, 
the  conductance  between  two  cubic  electrolyte  blocks  whose 
length  is  same  as  the  diameter  of  the  spherical  particle  ds  is 
necessary.  That  is 

°io-io  =  Kiods-  (8) 


3.  Calculation 

Geometrical  and  electrochemical  parameters  used  for  the 
random  resistor  network  simulations  are  summarized  in 
Table  1.  The  electrochemical  parameters  were  from  Sunde 
[2,3]  as  they  were  proven  to  produce  results  consistent  with 
published  experimental  observations.  For  example,  the  ionic 
conductivity  of  YSZ  and  the  electronic  conductivity  of  Ni 
are  standard  electrical  conductivities  measured  at  the  tem¬ 
perature  of  1000  °C  (0.1  S  cm"1  for  YSZ  and  2  x  104  S  cm'1 
for  Ni).  Only  the  composite  anode  composed  of  equal-sized 
ionic  and  electronic  particles  was  considered  due  to  the  limi¬ 
tation  of  the  regular  lattice.  However,  such  calculation  based 
on  the  regular  lattice  is  believed  to  be  useful  for  detecting 
important  features  of  the  transport  and  reaction  in  random 
particle  systems. 

The  numerical  simulations  were  conducted  by  solving  the 
current  conservation  equation  of  Eq.  (4)  for  each  percolated 
ionic  and  electronic  particles.  The  linearity  of  the  resistor  net¬ 
work  system  enables  the  arbitrary  choice  of  the  potentials  as 
boundary  conditions,  the  potentials  at  the  top  (pt op  and  at  the 
bottom  boundaries  <^bot-  By  prescribing  the  potential  differ- 


Table  1 


Geometric  and  electrochemical  parameters  for  random  resistor  network 
simulation 


Parameter 

Value 

Description 

ds 

2  fxm 

Diameter  of  spherical  particle 

1st  ly 

40  [xm 

Length  of  electrode  in  x-  and  ^-direction 

4 

50  |xm 

Thickness  of  electrode 

4 

20  [xm 

Thickness  of  electrolyte  layer 

4?  4> 

20  |xm 

Lengths  for  definition  of  penetrating 
electrolyte  structures  in  x-  and  y-direction 

4 

20  |xm 

Height  of  penetrating  electrolyte  structures 

hpb 

3  |xm 

TPB  length  per  contact,  /tpb  =  3  x  (0.5c/s) 

ho-iot  4l-el 

3  |xm 

Neck  perimeter  between  ions  and  electrons 

Kel 

2  x  104Scm‘’ 

Conductivity  of  electronic  phase  (Ni) 

^io 

0.1  S  cm1 

Conductivity  of  ionic  phase  (YSZ) 

/Cp 

10-4Scm‘1 

Polarization  conductivity  with  respect  to 
TPB  length 

0el 

0-1 

Volume  fraction  of  electronic  phase 

— yy  <pj)  <  (9) 

j  j 


where  the  convergence  criterion  s  was  10" 14  A. 

Once  the  potential  distribution  was  obtained,  the  currents 
flowing  through  the  top  and  the  bottom  boundaries  were  cal¬ 
culated  by  adding  currents  through  all  particles  adjacent  to 
the  boundaries.  The  total  current  flowing  through  the  com¬ 
posite  anode  /t  was  determined  by  averaging  the  currents 
calculated  at  the  top  and  the  bottom  boundaries.  The  total 
electrical  resistance  Rt  was  estimated  as 


Rt  T-  Rq 


^top  ^bot 

h 


1 

V 


(10) 


where  R&  is  an  electrical  resistance  of  an  electrolyte  layer 
which  was  included  to  reduce  the  effect  of  an  equi-potential 
condition  at  the  bottom  boundary.  In  the  simulation,  the  thick¬ 
ness  of  the  electrolyte  layer  was  fixed  to  20  pirn  (/e  =  20  p,m), 
and  Rq  was  calculated  by 


k 

k’io  Wy 


(ID 


Sunde  [2,3]  has  proposed  a  way  to  single  out  the  polariza¬ 
tion  resistance  Rv  from  the  total  resistance  Rt.  By  setting  the 
polarization  conductivity  to  be  infinite  (/cp  =  oo),  the  ohmic  re¬ 
sistance  Rq  was  evaluated  at  the  first  hand.  Then,  the  polariza¬ 
tion  resistance  Rv  was  determined  from  Rp  =  Rt  —  Rq,  where 
Rt  was  the  total  resistance  obtained  from  simulations  with 
an  appropriate  polarization  conductivity  /cp.  However,  it  is 
rather  insignificant  to  distinguish  the  polarization  resistance 
from  the  total  resistance  once  percolation  in  the  electronic 
particles  is  achieved.  When  the  volume  fraction  of  electronic 
phase  in  the  electrode  is  sufficiently  large  (0el  >  0-3),  the 
ohmic  resistance  Rq  becomes  small  compared  with  the  total 
resistance. 


4.  Results  and  discussions 

4.1.  Performance  enhancement  hy  penetrating 
electrolyte  structures 

Fig.  3  shows  the  calculated  area  specific  resistances  RVA  of 
the  composite  anodes,  with  and  without  the  penetrating  elec¬ 
trolyte  structures,  as  functions  of  the  volume  fraction  of  elec¬ 
tronic  phase  0el-  The  results  presented  as  symbols  in  Fig.  3 
are  averages  of  10  results  for  each  simulated  0el  with  dif¬ 
ferent  random  seed  numbers.  The  area  specific  polarization 
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Fig.  3.  Area  specific  polarization  resistances  RpA  (a)  and  total  resistances  RtA  (b),  with  respect  to  the  volume  fraction  of  electronic  phase  0ei-  Symbols  with 
solid  lines  represent  the  flat  electrolyte  (□),  and  the  electrolytes  with  penetrating  block  structures  (O),  rib  structures  (A),  and  cross-linked  rib  structures  (V). 
The  result  for  the  flat  electrolyte  by  Sunde  [2]  is  also  plotted  for  comparison  by  a  symbol  (■)  with  no  line. 


resistance  RpA  without  the  penetrating  electrolyte  structures, 
shown  in  Fig.  3(a),  is  in  good  agreement  with  the  result  by 
Sunde  [2] ,  where  small  discrepancy  can  be  partially  attributed 
to  the  dispersion  of  data  in  Monte-Carlo  simulation.  In  ad¬ 
dition,  the  present  calculations  excluded  the  particles  that  do 
not  belong  to  the  percolated  clusters,  and  this  approach  is 
different  from  that  of  Sunde  [2,3]. 

All  polarization  resistances,  shown  in  Fig.  3(a),  exhibit 
similar  trends  with  respect  to  0e l-  Starting  from  high  value, 
each  polarization  resistance  reaches  its  minimum  at  0.4  <  0ei 
<  0.5,  and  then  increases  again.  Such  behaviors  are  closely 
related  to  the  total  number  of  TPBs  XXpb  shown  in  Fig.  4(a). 
As  polarization  conductance  /cp  is  relatively  small  compared 
with  the  electronic  and  the  ionic  conductivities,  the  overall 
anodic  performance  is  governed  most  by  total  TPB  length 
X^tpb  (=XXpb  x  4pb)-  Therefore,  minimum  resistance  or 
maximum  performance  is  generally  observed  at  0el  ~  0.5 
where  X^tpb  becomes  largest.  Slight  discrepancy  between  0e i 
for  minimum  resistance  and  that  for  maximum  X^tpb  is  due  to 
large  difference  in  conductivities  of  the  ionic  and  electronic 
phases. 


Fig.  3(a)  shows  that  the  penetrating  electrolyte  structures 
reduce  the  polarization  resistance  and  thus  enhance  the  an¬ 
odic  performance,  especially  when  0ei  is  large.  As  noted  by 
Herbistritt  et  al.  [9],  one  reason  for  the  performance  enhance¬ 
ment  is  related  to  the  surface  enlargement  by  the  penetrating 
electrolyte  structures.  All  three  penetrating  electrolyte  struc¬ 
tures,  shown  in  Fig.  l(a)-(c),  were  found  to  double  the  sur¬ 
face  area  of  the  electrolyte  when  their  dimensions  were  set 
according  to  Table  1.  In  case  0el  >  0.7,  the  ionic  phase  cannot 
achieve  a  global  percolation  and  TPBs  are  spatially  restricted 
in  the  vicinity  of  the  electrolyte  surface.  That  results  in  XXpb 
directly  proportional  to  the  surface  area  of  the  electrolyte,  as 
shown  in  Fig.  4(a).  As  0el  approaches  1,  three  polarization 
resistance  curves  with  the  penetrating  electrolyte  structures 
fall  into  one  line,  about  half  of  that  with  the  flat  electrolyte 
(without  the  penetrating  electrolyte  structures). 

The  area  specific  total  resistances  RtA  are  shown  in 
Fig.  3(b)  for  0.35  <  0el  <  0.65,  where  percolation  of  both 
ionic  and  electronic  phases  is  ensured.  Minimum  RtA  is  ob¬ 
served  at  0.4  <  0ei  <  0.45  and  the  reduction  is  estimated  as 
much  as  20%  in  case  the  penetrating  block  structures  are 
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Fig.  4.  Total  number  of  TPBs  ^ntpb  with  respect  to  the  volume  fraction  of  electronic  phase  </>ei  (a),  and  the  spatial  distribution  of  TPBs  ntp b  for  <pe\  =  0.5  with 
respect  to  the  ion  path  length  Zp  (b).  Symbols  represent  the  flat  electrolyte  (□),  and  the  electrolytes  with  penetrating  block  structures  (O).  rib  structures  (A), 
and  cross-linked  rib  structures  (V). 
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considered.  Though  the  surface  enlargements  by  the  three 
penetrating  electrolyte  structures  are  same,  the  sacrifice  of 
the  total  TPB  length  X^tpb  due  to  the  addition  of  bulk  ionic 
phase  is  different,  which  causes  difference  in  the  anodic  per¬ 
formance  enhancement. 

For  a  closer  inspection  of  the  performance  enhancement 
by  the  penetrating  electrolyte  structures,  the  spatial  distribu¬ 
tion  of  TPBs  fttpb  as  a  function  of  the  relative  distance  from 
the  electrolyte  was  investigated.  For  this  purpose,  each  per¬ 
colated  ionic  particle  was  labeled  by  an  ion  path  length  /p 
which  was  defined  as  the  minimum  number  of  ionic  particles 
for  ions  to  travel  through  to  reach  the  bulk  electrolyte  phase, 
counting  itself  as  1.  Then,  the  number  of  TPBs  was  added 
for  the  ionic  particles  with  the  same  /p,  to  obtain  the  distri¬ 
bution  of  TPBs  ^tpb(^p)-  The  TPBs  formed  by  direct  contacts 
of  the  percolated  electronic  particles  onto  the  surface  of  the 
electrolyte  were  assigned  to  ntpb(0). 

Fig.  4(b)  suggests  that  the  penetrating  electrolyte  struc¬ 
tures  change  the  spatial  distribution  of  TPBs,  i.e.,  the  shape  of 
fttpb(/p);  they  increase  TPBs  that  have  shorter  ion  path  lengths 
at  the  sacrifice  of  TPBs  that  have  longer  ion  path  lengths. 
The  change  of  the  spatial  distribution  of  TPBs  is  another 
reason  for  the  performance  enhancement  by  the  penetrating 
electrolyte  structures.  The  three  penetrating  electrolyte  struc¬ 
tures  result  in  different  degrees  of  change;  largest  for  the 
penetrating  block  structures,  and  smallest  for  the  penetrating 
cross-linked  rib  structures.  Therefore  the  total  area  specific 
resistances  with  the  penetrating  block  structures  are  smallest 
as  shown  in  Fig.  3(b). 

The  current  generation  with  respect  to  the  ion  path  length 
/g(/p)  is  shown  in  Fig.  5(a).  It  is  observed  that  about  75%  of  the 
total  current  is  produced  by  TPBs  whose  ion  path  lengths  are 
less  than  10  particles  (/p  <  10)  in  case  of  the  flat  electrolyte. 
When  the  penetrating  electrolyte  structures  are  introduced, 
the  percentage  goes  up  to  85%.  By  dividing  the  current  gen¬ 
eration  by  the  number  of  TPBs  as  /g(/p)/«tpb(^p)  =  Apb(^p)* 
the  current  generation  efficiency  of  TPBs  as  a  function  of  the 
ion  path  length  is  investigated.  As  shown  in  Fig.  5(b),  the 
current  generation  efficiency  /tpb(4>)  is  an  exponentially  de- 


4  (fim) 


Fig.  6.  Area  specific  polarization  resistances  RpA  with  respect  to  the  elec¬ 
trode  thickness  lz  for  </>ei  =  0-5.  Symbols  represent  the  flat  electrolyte  (□), 
and  the  electrolytes  with  penetrating  block  structures  (OX  rib  structures  (A), 
and  cross-linked  rib  structures  (V). 

caying  function  that  decreases  about  50%  at  every  increase  of 
the  ion  path  length  by  5  particle  unit.  That  is  to  say,  TPBs  at 
the  electrolyte  surface  produce  about  four  times  more  cur¬ 
rent  than  TPBs  of  the  ion  path  length  of  10  particle  unit 
(4pb(0)  ~  4/tpb(10)).  As  TPBs  that  have  shortest  ion  path 
length  are  most  effective,  rather  small  difference  in  perfor¬ 
mance  enhancement  by  three  penetrating  electrolyte  struc¬ 
tures,  shown  in  Fig.  3(b),  can  be  explained. 

Variation  of  the  area  specific  polarization  resistance  RpA 
according  to  the  thickness  of  the  composite  electrode  lz  at  0ei 
=  0.5  is  shown  in  Fig.  6.  The  calculated  polarization  resis¬ 
tance  for  the  flat  electrolyte  shows  a  good  agreement  with  the 
theoretical  prediction  by  Kenjo  et  al  [1 1],  which  is  expressed 
as 

Rp=(f  =  fpkcoth<ff  (12) 

1  o  V  k/ 

where  (po  and  7o  are  the  potential  difference  and  the  current,  p 
is  the  resistivity  of  the  electrolyte,  and  /  is  the  electrode  thick¬ 
ness.  In  Fig.  6,  all  polarization  resistances  initially  decrease 


Fig.  5.  Distribution  of  current  generation  7g  (a)  and  the  current  generation  per  TPB  7g/ntpb  (b),  for  0e l  =  0.5  with  respect  to  the  ion  path  length  /p.  Symbols 
represent  the  flat  electrolyte  (□),  and  the  electrolytes  with  penetrating  block  structures  (OX  rib  structures  (A),  and  cross-linked  rib  structures  (V). 
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Fig.  7.  Area  specific  polarization  resistances  RpA  (a)  and  change  of  the  spatial  distribution  of  TPBs  (b),  with  respect  to  the  heights  of  the  penetrating  electrolyte 
structures  /c,  for  0ei  =  0.5  and  lz  =  50  fxm.  Symbols  represent  the  electrolytes  with  penetrating  block  structures  (OX  rib  structures  (A),  and  cross-linked  rib 
structures  (V). 


with  the  increase  of  the  anode  thickness  and  then  remain 
rather  constant.  When  the  thickness  of  the  composite  anode 
is  small,  TPBs  in  the  anode  are  limited  in  number;  this  results 
in  large  polarization  loss  and  resistance.  However,  once  suffi¬ 
cient  number  of  TPBs  is  ensured,  an  additional  increase  of  the 
anode  thickness  does  not  reduce  the  polarization  resistance. 
For  the  flat  electrolyte,  TPBs  in  the  anode  become  sufficient 
for  the  electrochemical  reaction  when  the  anode  thickness  lz 
is  larger  than  20  p,m.  In  Fig.  5(a),  the  current  produced  by 
TPBs  that  have  less  than  20  particle  unit  (Zp  <  20)  amounts 
more  than  90%  of  total  current. 

4.2.  Optimum  design  of  penetrating  electrolyte 
structures 

As  a  first  step  for  the  optimization  of  the  penetrating  elec¬ 
trolyte  structures,  the  behaviors  of  the  area  specific  polar¬ 
ization  resistances  RpA  of  the  composite  anode  with  varying 
height  lc  are  shown  in  Fig.  7(a).  The  anodic  polarization  re¬ 


sistance  with  the  penetrating  rib  structures  reaches  minimum 
at  lc  ~  40  p,m  while  that  with  the  penetrating  block  structures 
remains  constant  for  lc  >  40  jim.  The  polarization  resistance 
with  the  penetrating  cross-linked  rib  structures  shows  a  dis¬ 
tinctive  minimum  at  lc  ~  35  pm. 

The  distributions  of  TPBs  with  respect  to  the  ion  path 
length  lp  are  shown  in  Fig.  7(b),  for  different  heights  lc 
of  the  penetrating  block  structures.  Increasing  the  height  of 
the  penetrating  block  structures  continuously  increases  TPBs 
whose  ion  path  length  is  smaller  than  10  particle  unit  while 
decreasing  TPBs  that  have  the  longest  ion  path  length.  Since 
TPBs  with  smaller  ion  path  length  are  more  efficient  in  elec¬ 
trochemical  reaction  and  current  generation,  the  anodic  per¬ 
formance  increases  with  the  increase  of  the  height  of  the 
penetrating  block  structures.  The  spatial  distribution  of  TPBs 
shows  that  the  ion  path  lengths  of  all  TPBs  are  within  20  par¬ 
ticle  unit  for  the  block  height  of  32  pm,  and  within  10  particle 
unit  for  the  block  height  of  40  pm.  Then,  further  increase  of 
the  height  of  the  penetrating  block  structures  cannot  improve 


Fig.  8.  Area  specific  polarization  resistance  RpA  of  the  composite  anode  with  penetrating  block  structures  for  $ei  =  0-5.  The  heights  of  the  penetrating  block 
structures  lc  are  30  |xm  (a)  and  40  |xm  (b),  and  the  thickness  of  composite  anode  lz  is  50  |xm. 


28 


D.H.  Jeon  et  al.  /  Journal  of  Power  Sources  139  (2005)  21-29 


Fig.  9.  Area  specific  polarization  resistance  RpA  of  the  composite  anode  with  penetrating  rib  structures  for  0ei  =  0-5.  The  heights  of  the  penetrating  rib  structures 
/c  are  30  |xm  (a)  and  40  |xm  (b),  and  the  thickness  of  composite  anode  lz  is  50  |xm. 


Fig.  10.  Area  specific  polarization  resistance  RpA  of  the  composite  anode  with  penetrating  cross-linked  rib  structures  for  0ei  =  0-5.  The  heights  of  the  penetrating 
cross-linked  rib  structures  /c  are  30  |xm  (a)  and  40  |xm  (b),  and  the  thickness  of  composite  anode  lz  is  50  |xm. 


the  anodic  performance  and  may  deteriorate  the  performance 
on  the  contrary. 

From  Fig.  7,  the  criterion  for  the  height  of  penetrating 
electrolyte  structures  (for  50  p,m  thick  composite  anode)  is 
identified  and  the  heights  of  30  and  40  p,m  are  chosen  for  fur¬ 
ther  optimization.  The  polarization  resistances  obtained  by 
varying  the  planar  dimensions  of  the  penetrating  electrolyte 
structures,  lc  and  4>,  are  shown  in  Figs.  8-10.  It  is  found  that 
about  40-45%  reduction  in  the  anodic  polarization  resistance 
can  be  achieved  by  the  optimization  of  the  planar  dimensions 
and  the  difference  in  the  minimum  resistances  with  three  pen¬ 
etrating  electrolyte  structures  is  rather  small;  99  m  Q  cm2  for 
the  penetrating  block  and  rib  structures,  and  95  m  Q  cm2  for 
the  penetrating  cross-linked  rib  structures,  compared  with 
184  m  Q  cm2  for  the  flat  electrolyte.  The  minimum  polariza¬ 
tion  resistances  are  also  rather  insensitive  to  the  height  of 
the  penetrating  electrolyte  structures,  resulting  only  about 
5%  difference  for  the  height  of  30  and  40  p,m.  The  optimum 


anodic  performance  with  the  penetrating  block  structures  is 
achieved  when  /a  >  4>,  while  that  with  the  penetrating  cross- 
linked  rib  structures  is  achieved  when  /a  <  4>.  The  penetrating 
rib  structures  are  found  to  perform  best  when  /a  ~  Zb-  When 
the  ratio  4,//a  is  equal,  the  penetrating  electrolyte  structures 
with  smaller  /a  +  4>  perform  better,  for  the  surface  enlarge¬ 
ment  is  larger. 

5.  Conclusions 

Attachment  of  the  penetrating  electrolyte  structures  to  the 
conventional  flat  electrolytes  of  SOFCs  was  proposed,  and  the 
reduction  in  the  anodic  polarization  resistances  was  quantita¬ 
tively  assessed  using  a  random  resistor  network  model.  The 
predicted  dependency  of  the  anodic  performance  on  micro- 
structural  parameters  is  consistent  with  the  results  in  the 
literature,  showing  that  the  anodic  performance  is  mainly 
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governed  by  the  reaction  efficiency  and  the  total  TPB  length 
of  the  composite  anode.  The  best  anodic  performances  were 
obtained  when  the  volume  fraction  of  electronic  phase  is  be¬ 
tween  0.4  and  0.5,  and  when  the  thickness  of  the  anode  is 
sufficiently  large. 

The  simulation  results  demonstrated  that  the  penetrating 
electrolyte  structures  could  reduce  the  anodic  polarization 
resistance  considerably.  From  the  investigation  on  the  dis¬ 
tribution  of  TPBs,  the  shift  in  the  spatial  distribution  due 
to  the  penetrating  electrolyte  structures  is  identified  as  the 
main  mechanism  for  the  anodic  performance  enhancement. 
While  total  TPB  length  is  reduced  by  the  introduction  of 
the  penetrating  electrolyte  structures,  TPBs  that  are  close  to 
the  electrolyte  surface  are  increased  at  the  sacrifice  of  TPBs 
that  are  farthest  from  the  electrolyte  surface.  As  the  current 
generation  efficiency  of  TPBs  is  found  to  be  the  exponen¬ 
tially  decreasing  function  of  the  distance  from  the  electrolyte 
surface,  such  change  in  the  spatial  distribution  of  TPBs  is 
favorable  for  the  anodic  performance. 

A  parametric  study  was  conducted  to  derive  design  crite¬ 
ria  for  the  dimensions  of  the  penetrating  electrolyte  structures 
to  obtain  optimal  anodic  performance.  Larger  height  of  the 
penetrating  electrolyte  structures  is  favorable  for  the  anodic 
performance  by  facilitating  easy  penetration  of  oxygen  ions 
into  the  composite  anode  as  well  as  by  increasing  the  surface 


area  of  the  electrolyte  more.  Optimal  heights  of  the  penetrat¬ 
ing  electrolyte  structures  are  estimated  to  be  between  60  and 
80%  of  the  thickness  of  the  composite  anode.  As  to  the  opti¬ 
mal  ranges  for  the  planar  dimensions,  minimum  polarization 
resistances  are  achieved,  when  /a  >  4>  for  the  penetrating  block 
structures,  when  /a  <  4>  for  the  penetrating  cross-linked  rib 
structures,  and  when  /a  ~  4>  for  the  penetrating  rib  structures. 

References 

[1]  N.Q.  Minh,  J.  Am.  Ceram.  Soc.  76  (1993)  563. 

[2]  S.  Sunde,  J.  Electrochem.  Soc.  142  (1995)  L50. 

[3]  S.  Sunde,  J.  Electrochem.  Soc.  143  (1996)  1930. 

[4]  J.  Abel,  A.A.  Kornyshev,  W.  Lehnert,  J.  Electrochem.  Soc.  144 
(1997)  4253. 

[5]  P.  Costamagna,  P.  Costa,  V.  Antonucci,  Electrochim.  Acta  43  (1998) 
375. 

[6]  P.  Costamagna,  M.  Panizza,  G.  Cerisola,  A.  Barbucci,  Electrochim. 
Acta  47  (2002)  1079. 

[7]  S.H.  Chan,  Z.T.  Xia,  J.  Electrochem.  Soc.  148  (2001)  A388. 

[8]  X.J.  Chen,  S.H.  Chan,  K.A.  Khor,  Electrochim.  Acta  49  (2004)  1851. 

[9]  D.  Herbistritt,  A.  Weber,  E.  Ivers-Tiffee,  J.  Eur.  Ceram.  Soc.  21 
(2001)  1813. 

[10]  V.T.  Srikar,  K.T.  Turner,  T.Y.  Andrew  Ie,  S.M.  Spearing,  J.  Power 
Sources  125  (2004)  62. 

[11]  T.  Kenjo,  S.  Osawa,  K.  Fujikawa,  J.  Electrochem.  Soc.  138  (1991) 
349. 


